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METHOD FOR IDENTIFICATION OF AN OSCILLATION IN AN 
ELECTRICAL POWER SUPPLY SYSTEM 

CLAIM FOR PRIORITY 
5 This application claims priority to International 

Application No. PCT/DEOO/03434 which was filed in the 
German language on September 11, 2000. 

TECHNICAL FIELD OF THE INVENTION 
10 The invention relates to a method for producing at 

least one signal (oscillation signal) , which indicates 
an oscillation in an electrical power supply system. 

BACKGROUND OF THE INVENTION 
15 German Laid-Open Specification DE 195 03 626 Al 

describes a method of identifying an oscillation. In 
this method, once a memory element has been set, 
further impedance values are checked to determine 
whether the oscillation that has been found is still 
20 continuing, by determining the. rate of change of the 
magnitude of respectively successive impedance values. 
If it is found that the rate of change is above a limit 
value, this identifies that the oscillation has 
stopped, and the memory element is reset. It is 
25 difficult to define such a limit value, particularly 
when a large number of generators are connected in the 
power supply networks, and * complex oscillations can; 
thus occur. 

SUMMARY OF THE INVENTION 
30 The invention relates to a method for producing at 

least one signal (oscillation signal) , which indicates 
an oscillation in an electrical power supply system, in 
which method the phase current and the phase voltage 
are in each case sampled from at least one phase of the 
35 power supply system, forming phase current and phase 
voltage sample values, impedance values are formed from 
the phase current and phase voltage sample values, the 
impedance values are monitored for the presence of any 
oscillation and, if an oscillation is identified, at 



least one memory element is set, and the oscillation 
signal is output at its output, after setting the 
memory element, further impedance values are checked to 
determine whether the oscillation that has been found 
5 is still continuing, the. memory element remains 
uninfluenced if the oscillation continues, and the 
memory element is reset if the oscillation has stopped. 

The invention specifies a method of detecting the 
oscillation behavior of an electrical power supply 

10 system in a safe and reliable manner. 

In one embodiment of the invention, there is a 
method that checks the impedance values and makes use 
of an oscillation model which is formed from previous 
impedance values associated with the oscillation, or 

15 from variables which are dependent on these impedance 
values. A check is then carried out to determine 
whether a further impedance value formed at. that time 
or a variable which is dependent on the further 
impedance value differs from the oscillation model, and 

20 any occurrence of a further impedance value or of a 
variable dependent on this impedance value which 
differs from the oscillation model is assessed as the 
oscillation having stopped. 

One advantage of the method according to the 

25 invention is that the oscillation model allows even 
complex oscillations to be described, and it is thus 
possible to identify that the oscillation has stopped 
with a high level of reliability even in the case of 
such complex oscillations. 

30 The oscillation model can advantageously be 

determined by means of a least squares estimation 
method. This estimation method allows a mathematical 
oscillation model to be produced from successive 
impedance values which have been formed after the 

35 setting of the memory element, that is to say after the 
start of the oscillation. 

A function in the form f (x) = ax^+bx^+cx+d with the 
parameters a, b, c and d can be used as the model rule 
for the oscillation model, in which one or more 
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parameters can be defined to be zero before the start 
of the estimation method. First, second or third order 
power functions can be used as the model rule. 
Furthermore, a sum of sine and cosine functions, which 
5 decay with time, can be used as the model rule for the 
oscillation model. These model rules make it possible 
to describe even complex oscillations mathematically. 

The oscillation model can be formed directly for 
the determined impedance values of the oscillation, or 

10 else for variables dependent on these impedance values. 
Resistance values. R, reactance values X, time 
derivative values dZ/dt of the impedance, time 
derivative values dR/dt of a resistance or time 
derivative values dX/dt of a reactance can be used as 

15 dependent variables. Choice of the most suitable 
variable for the oscillation model makes it possible to 
determine with a high level of reliability that the 
oscillation has stopped, with the choice of the 
variable being dependent on the individual system 

20 configuration of the electrical power supply system. 

In "one advantageous embodiment of the invention, 
positive phase sequence system impedance values can be 
formed from the phase current and phase voltage sample 
values, and a common memory element can be provided, 

25 and a common oscillation signal can be produced, for 
the phases in the power supply system. This variant can 
be used when the aim is to make a statement relating to 
any oscillation occurring at the same time in all the 
phases in the power supply system. 

30 In a further embodiment of the method according to 

the invention, phase impedance values are formed from 
the phase current and phase voltage sample values of 
each phase of the power supply system to be 
investigated for oscillation, and a dedicated memory 

35 element is provided, and a dedicated oscillation signal 
is produced, for each of these phases. In this 
embodiment, the oscillation response of each individual 
phase in the power supply system can be investigated 
separately. That is, both the starting and the 
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stopping of an oscillation are identified. This is 
particularly advantageous when oscillations occur in a 
single phase, but not in all the phases, in the power 
supply system. Oscillations such as these frequently 
5 occur in the case of so-called single-pole pauses in 
high-voltage systems. Single-pole pauses are produced 
by single-pole conductor ground faults, which can be 
expected frequently in high-voltage systems, and in 
which an arc is struck between one conductor and 
10 ground. In this type of fault, a single-pole pause is 
produced. That is, the single phase in which the 
single-pole conductor -ground fault has occurred is 
switched off briefly. The arc is thus quenched, and the 
fault is frequently corrected. Switching off a single 
15 pole of one phase can result in oscillations occurring 
in the remaining phases which are not switched off. 
These oscillations cannot be identified, for example, 
by monitoring the positive phase sequence system 
impedance values, since positive phase sequence system 
20 impedance values are formed when sample values are 
available for the phases in the power supply system- In 
the case of a single-pole pause, it is advantageous to 
be able to produce a dedicated oscillation signal for 
each phase in the power supply system. During the 
25 single-pole pause, this oscillation signal is produced 
for those phases which are not switched off. The 
oscillation behavior of the power supply system can 
thus be determined individually for each phase, and 
independently of the state of the other phases. 
3 0 The phase impedance values of the individual 

phases in the electrical power supply system can, for 
example, be formed by, in order to form the phase 
impedance values, 

a variable U_re including the real part of the 
35 phase voltage sample values, a variable U_im including 
the imaginary part of the phase voltage sample values, 
a variable I_re including the real part of the phase 
current sample values and a variable I_im including the 
imaginary part of the phase current sample values being 



formed from the phase current and phase voltage sample 
values (i, u) for each phase, 

a phase real power variable P being determined 
from P = U_re-I_re - U_iTn-I_im, 
5 - a phase wattless component variable Q being 
determined from Q = U_im-I_re + U_im-I_re, 
a squared phase current amplitude variable 
being determined from = I_re-I_re + I_im-I_im, 
system- frequency components in each case being 
10 removed by means of least squares estimation 

method from the phase real power variable P, from 
the phase wattless component variable Q and from 
the squared phase current amplitude variable 1^ , 
and 

15 - . phase resistance values R being determined from 
R=P/I^ and phase reactance values X being 
determined from X=Q/I^, and phase impedance values 
Z=R+jX being thus determined. 

20 When forming the phase impedance values, it is 

advantageous to remove the system frequency components 
(for example 50 Hz components) from the phase real 
power variable from the phase wattless component 

variable Q and from the squared phase current amplitude 

25 variable by means a least squares estimation method 
in each case. Such system frequency interference 
components would adversely affect the evaluation of the 
phase impedance values determined from these variables. 

30 BRIEF DESCRIPTION OF THE DRAWINGS 

In order to explain the invention further. 
Figure 1 shows a block diagram of an exemplary 
embodiment of the method according to the invention. 

Figure 2 shows a block diagram for determining the 
35 phase impedance values. 

Figure 3 shows the impulse response of the filters 
used for impedance determination. 

Figure 4 shows the real power and wattless 
component variable profiles before filtering. 



Figure 5 shows the real power and wattless 
component variable profiles after filtering. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
Figure 1 shows, schematically, a method for 
determining the oscillation behavior of a three-phase 
electrical power supply, system, by means of which a 
dedicated oscillation signal Pdl , Pd2 and Pd3 is 
produced for each phase in the power supply system. 
This is done by providing three changeover switches Ul, 
U2 and U3 and three memory elements Spi, Sp2 and Sp3 . 
The connecting lines between the individual units in 
the layout are designed with three poles. The phase 
current and phase voltage sample values i and u of all 
three phases are supplied to a unit for impedance 
determination lb, at whose output phase impedance 
values Z are output for the three phases. These phase 
impedance values Z are supplied via the changeover 
switches Ul, U2 and U3 to an oscillation identification 
unit Pe. The oscillation identification unit Pe uses 
the time profiles of the phase impedance values Z to 
identify the occurrence of any oscillation in the 
individual phases, for example in phase 1, and emits an 
oscillation set signal Ps at its output for each phase 
in which oscillation is. identified, for example for the 
phase 1. The oscillation set signal Ps sets the memory 
element associated with the respective phase, for 
example Spl, and this memory element emits at its 
output the phase-specific oscillation signal, for 
example, Pdl . When an oscillation signal is being 
emitted, for example the oscillation signal Pdl, the 
changeover switch, for example Ul, associated with the 
respective phase is switched over. The phase impedance 
values Z which are still formed for the phase in which 
the oscillation has been identified, for example the 
phase 1, are supplied to an oscillation signal 
resetting unit Pu. This oscillation signal resetting 
unit Pu identifies that the oscillation has stopped 
and, in this case, emits an oscillation reset signal Pr 



at its output, which resets the memory element for the 
respective phase, for example Spl , The oscillation 
signal for the respective phase, for example Pdl,. is 
thus also no longer emitted, and the respective 
changeover switch, for example Ul, moves back to its 
original position once again. In response to a 
stimulus, a unit for phase selection Pa ensures that 
the phase impedance values of the phases to be 
investigated for oscillation are in each case processed 
by the oscillation identification unit Pe and by the 
oscillation signal resetting unit Pu. 

The method of operation of the four units 
comprising impedance determination lb, the oscillation 
identification unit. Pe, the oscillation signal 
resetting unit Pu and phase selection Pa will be 
explained in more detail in the following text. 

As shown in Figure 2, the phase current and phase 
voltage sample values i and u are filtered in the 
impedance determination unit lb by means of orthogonal 
FIR filters Fl, F2, F3 and F4, thus resulting in the 
production of a variable U_re containing the real part 
of the phase voltage sample values, a variable U_im 
containing the imaginary part of the phase voltage 
sample values, a variable I_re containing the real part 
of the phase current sample values, and the variable 
I_im containing the imaginary part of the phase current 
sample values. 

Figure 3 shows the impulse responses of the 
filters Fl to F4, with the impulse response of the 
filters Fl and F3 which determine the real parts being 
annotated "o", and the impulse response of the filters 
F2 and F4 which determine the imaginary parts, being 
annotated . 

As shown in Figure 2, following this, a phase real 
power variable P is calculated in accordance with 
equation (1) below and a phase wattless component 
variable Q is calculated in accordance with equation 
(2) in the unit 5, and a squared phase current 
amplitude variable is calculated in accordance, with 



equation (3) in the unit 6. 

p = U_re-I_re - U_im-I_um (1) 

Q = U_ini-I_re + U_im'I__re (2) 

= I re-I re + I im-I im (3) 



After this, the phase real power variable P, the 
5 phase wattless component variable Q and the scfuared 
phase current amplitude variable are filtered in 

units 7 and 8 in order to remove the interference 50 Hz 
components contained in these variables; this results 
in the filtered variables P' , Q' and I^' . The least 
10 squares estimation method used for this filtering will 
be explained in detail further below. 

The upper illustration a) in Figure 4 shows the 
profile of the real power variable P, and the lower 
illustration b) shows the profile of the Wattless 
15 component variable Q, before filtering by means of the 
least squares estimation method, in each case plotted 
against the time, t . 

The upper illustration a) in Figure 5 shows the 
profile of the real power variable P' , and the lower 
20 illustration b) shows the profile, of the Wattless 
component variable Q' , after filtering by means of the 
least squares estimation method. It can clearly be 
seen that the 50 Hz components have been removed. 

As shown in Figure 2, after the filtering in the 
25 unit 9, phase resistance values R and phase reactance 
values X are determined in accordance with equation 
(4) , and the phase impedance values Z=R+jX determined 
from them are emitted at the output of the impedance 
determination unit lb. 

30 

R==P'/I^' X=:Q'/I^» (4) 

A least squares estimation method using a signal 
model in accordance with equation (5) is applied 
. separately to each of the variables P, Q and in 
35 order to filter out the 50 Hz components included in 
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the phase real power variable P, in the phase Wattless 
component variable Q and in the squared phase current 
amplitude variable I^- 

y, =A c ^•sin(o)ok T^) + B e * .cos(cOok T^) + C (S) 

5 

The estimation method uses the squared phase 
current amplitude variable I^, the phase real power 
variable P and the phase Wattless component variable Q 

10 to calculate the parameters A, B and C in the signal 
model. The parameter C gives the sought magnitude of 
the phase real power variable P', of the phase Wattless 
component variable Q' and of the squared phase current 
amplitude variable ' . The summands with the 

15 parameters A and B model the 50 Hz components. The 
variable (oo is the frequency (50 Hz) to be filtered 
out, and Ta is the sampling time. 

If an equivalent circuit with two generator 
machines at the ends of a power transmission line is 

20 considered for the power supply system, the amplitude 
of the 50 Hz components decays with the time constants 
X of the sum impedance between the two generator 
machines in accordance with equation (6) , where L is 
the loop inductance and R is the loop resistance of the 

25 circuit which is closed via the two generator machines. 

The coefficients A, B and C are determined such 
30 that the sum of the squares of the errors between 
values y determined from the phase current and phase 
voltage sample values i and u and the sample values yk 
calculated in accordance with equation (5) becomes a 
minimum (see equation (7)). 



35 



(7) 



In equation (7) , J represents the Q-criterion to 
be minimized. The signal model included in equation (5) 
5 is used as the function h(0k). The parameters A, B and 
C to be determined form a vector 0jc in accordance with 
equation (8) . 



10 



15 



B 



(8) 



The Q-criterion J is derived based on the 
parameter vector 0k in order to solve the minimization 
task. For the signal model in accordance with equation 
(5), this, then results in equation (9) together with 
equation (10) . 



0= S2x^(yi-Y ©0 

l=fc-N 



(9) 



* 

Y = 



cos^ 



(10) 



If equation (9) is solved for the vector &x, then 
this results in equation (11) , by means of which, and 
using the matrix Sj^ included in equations (12) and 
(13) , the vector ©jc is determined. 



Sk = 



(13) 



Of the parameters A, B and C. included in the 
vector ©kr only the parameter C is evaluated. The 
10 vectors in accordance with equation (10) and the 

matrix Sk in accordance with equation (13) are 
calculated and are stored as constants, so that they 
are available when the method is used. 

Monotony criteria are applied to the locus curves 

15 of the impedance values in the impedance plane in the 
oscillation identification unit Pe, in order to 
identify the oscillation process. This method for 
identification of the oscillation process is described 
in German Patent DE 197 46 719 CI. 

20 The oscillation signal resetting unit Pii 

determines whether an oscillation which has already 
been identified is still continuing. The procedure used 
for this purpose comprises the production of an 
oscillation model for phase impedance values Z 

25 associated with the oscillation. A check is then 
carried out to determine whether the locus curve which 
is described by the newly determined phase impedance 
values Z still corresponds to the oscillation model. 
When producing the oscillation model, it is assumed 



that the locus curve is free of discontinuities, and 
that its direction changes slowly. In the present 
exemplary embodiment, the locus curve is described by a 
first order power function, that is to say by a linear 
equation, in accordance with equation (14) . 

X(R)=m-R+Xo (14) 

The parameters m and XO are determined by means of 
a non-recursive least squares estimation method from 
the last N determined phase impedance values Z. 

The linear equation is used as a model rule for 
the least squares estimation method, with the parameter 
m characterizing the gradient, and the parameter XO 
characterizing the offset of the linear equation. The 
parameters m and XO for the model in accordance with 
equation (14) are determined from the last determined 
value pairs (Ri, Xi) of the phase impedance values Zi 
such that the sum of the squares of the errors between 
the values Xi determined from the measured phase 
current and phase voltage sample values i. and u and the 
values X calculated in accordance with equation (14) is 
a minimum (see equation (15)) . 

J' t(^,-h(eJ)'-»MIN (15) 

In equation (15) , J is the Q-criterion to. be 
minimized, the model rule in accordance with equation 
(14) is used as the function h(0k> . In accordance with 
equation (16) , the parameter vector 0^ contains the 
parameters m and XO, to be determined, from the model 
rule . 



In order to solve the minimization task, the Q- 
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criterion J is derived on the basis of the parameter 
vector 0K. The following equations (17) and (18) are 
then obtained for the signal model in accordance with 
equation (14) : 

0= i2r;(x.-/.(ej,) (17) 



(18) 



If equation (17) is solved for the parameter 

vector ©k, then this results in equation (19) for 

10 determining the parameter vector ©j,. 

where 



(21) 



After substitution of the parameters in equation 
(14), this results in the estimated oscillation model. 
If a newly determined phase impedance value Z 
corresponds to the oscillation model. That is, it is 

20 within a tolerance band around the linear equation 
represented by equation (14) , then this identifies that 
the oscillation is continuing. If the newly determined 
phase impedance value Z is outside the tolerance band, 
then this identifies that the oscillation has stopped, 

25 and an oscillation reset signal Pr is emitted at the 
output of the oscillation signal resetting unit Pu for 
the memory element Spl, Sp2 or Sp3 for the respective 
phase - 

The phase selection unit Pa receives a stimulus 
30 from, . for example, a distance protection device which 
is not illustrated. Depending on the nature of the 



stimulated loop, it determines the phases for which the 
oscillation identification unit Pe and/or the 
oscillation signal resetting unit Pu should investigate 
the oscillation behavior. The following table shows the 
association between the stimulated loops and the 
phases . 



Stimulated loops 


Phases to be investigated 
for oscillation behavior 


LIE 


LI 


L2E 


L2 


L3E 


L3 


L12 


LI and L2 


L23 


L2 and L3 


L31 


LI and L3 
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Description 

Method for identification of an oscillation in an 
electrical power supply system 

5 

The invention relates to a method for producing at 
least one signal (oscillation signal) , which indicates 
an oscillation in an electrical power supply system, in 
which method the phase current and the phase voltage 

10 are in each case sampled from at least one phase of the 
power supply system, forming phase current and phase 
voltage sample values, impedance values are formed from 
the phase current and phase voltage sample values, the 
impedance values are monitored for the presence of any 

15 oscillation and, if an oscillation is identified, at 
least one memory element is set, and the oscillation 
signal is output at its output, after setting the 
memory element, further impedance values are checked to 
determine whether the oscillation that has been found 

20 is still continuing, the memory element remains 
uninfluenced if the oscillation continues, and the 
memory element is reset if the oscillation has stopped . 

A method such as this is described in German Laid-Open 
25 Specification DE 195 03 626 Al . In this method, once 
the memory element has been set, further impedance 
values are checked to determine whether the oscillation 
that has been found is still continuing, by determining 
the rate of change of the magnitude of respectively 
30 successive impedance values and, if it is found that 
the rate of change is above a limit value, this 
identifies that the oscillation has stopped, and the 
memory element is reset. It has been found to be 
difficult to define such a limit value, particularly 
3 5 when a large number of generators are connected in the 
power supply networks, . and complex oscillations can 
thus occur. 
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The invention is based on the object of specifying a 
method by means of which the oscillation behavior of an 
electrical power supply system can be detected safely 
and reliably at all times. 

For a method of the type mentioned initially, this 
object is achieved according to the invention in that 
the check of the further impedance values makes use of 
an oscillation model which is formed from previous 
impedance values associated with the oscillation, or 
from variables which are dependent on these impedance 
values; a check is then carried out to determine 
whether a further impedance value formed at that time 
or a variable which is dependent on this further 
impedance value differs from the oscillation model, and 
any occurrence of a further impedance value or of a 
variable dependent on this impedance value which 
differs from the oscillation model is assessed as the 
oscillation having stopped. 

One major advantage of the method according to the 
invention is that the oscillation model allows even 
complex oscillations to be described, and it is thus 
possible to identify that the oscillation has stopped 
with a high level of reliability even in the case of 
such complex oscillations. 

The oscillation model can advantageously be determined 
by means, of a least squares estimation method. This 
estimation method allows a mathematical oscillation 
model to be produced from successive impedance values 
which have been formed after the setting of the memory 
element, that is to say after the start of the 
oscillation, 

A function in the form f (x) = ax^+bx^+cx+d with the 
parameters a, b, c and d can be used as the model rule 
for this oscillation model, in which one or more 
parameters can be defined to be zero before the start 
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of the estimation method. First, second or third order 
power functions can thus be used as the model rule. 
Furthermore, a sum of sine and cosine functions, which 
decay with time, can be used as the model rule for the 
5 oscillation model. These model rules make it possible 
to describe even complex oscillations mathematically* 

The oscillation model can be formed directly for the 
determined impedance values of the oscillation, or else 

10 for variables dependent on these impedance values. 
Resistance values R, reactance values X, time 
derivative values dZ/dt of the impedance, time 
derivative values dR/dt of a resistance or time 
derivative values dx/dt of a reactance can be used as 

15 dependent variables . Choice of the most suitable 
variable for the oscillation model makes it possible to 
determine with a high level of reliability that the 
oscillation has stopped, with the choice of the 
variable being dependent on the individual system 

20 configuration of the electrical power supply system. 

In one advantageous embodiment of the invention, 
positive phase sequence system impedance values can be 
formed from the phase current and phase voltage sample 

25 values, and a common memory element can be provided, 
and a common oscillation signal can be produced, for 
all the phases, in the power supply system. This variant 
can be used when the aim is to make a statement 
relating to any oscillation occurring, at the same time 

30 in all the phases in the power supply system. 

In a further advantageous embodiment of the method 
according to the invention, phase impedance values are 
formed from the phase current and phase voltage sample 
35 values of each phase of the power supply system to be 
investigated for oscillation, and a dedicated memory 
element is provided, and a dedicated oscillation signal 
is produced, for each of these, phases. In this 
embodiment, the oscillation response of each individual 



phase in the power supply system can be investigated 
separately, that is to say both the starting and the 
stopping of an oscillation are identified. This is 
particularly advantageous when oscillations occur in 
5 only a single phase, but not in all the phases, in the 
power supply system. Oscillations such as these 
frequently occur in the case of so-called single-pole 
pauses in high-voltage systems. Single-pole pauses are 
produced by single -pole conductor ground faults, which 
10 can be expected frequently in high-voltage systems, and 
in which an arc is struck between one conductor and 
ground. In this type of fault, a single-pole pause is 
produced that is to say the single phase in which the 
single-pole conductor-ground fault has occurred is 
15 switched off briefly. The arc is thus quenched, and the 
fault is frequently corrected. Switching off a single 
pole of one phase can result in oscillations occurring 
in the remaining phases which are not switched off. 
These oscillations cannot be identified, for example, 
20 by monitoring the positive phase sequence system 
impedance values, since positive phase sequence system 
impedance values can be formed only when sample values 
are available for all the phases in the power supply 
system. In the case of a single-pole pause, it is now 
25 highly advantageous to be able to produce a dedicated 
oscillation signal for each phase in the power supply 
system; during the single-pole pause, this oscillation 
signal is produced only for those phases which are not 
switched off. The oscillation behavior of the power 
30 supply system can thus be determined individually for 
each phase, and independently of the state of the other 
phases . 

The phase impedance values of the individual phases in 
35 the electrical power supply system can, for example, be 
formed by, in order to form the phase impedance values, 

a variable U_re containing the real part of the 
phase voltage sample values, a variable, U_im containing 
the imaginary part of the phase voltage sample values. 
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a variable I_re containing the real part of the phase 
current sample values and a variable I_im containing 
the imaginary part of the phase current sample values 
being formed from the phase current and phase voltage 
5 sample values (i, u) for each phase, 

a phase real power variable P being determined 
from P = U_re-I_re - U_im-I_im 

a phase wattless component variable Q being 
determined from Q = U_im-I_re + U__im-I__re 

10 - a squared phase current amplitude variable 
being determined from = I_re-I_re + I_im-I_im 
system- frequency components in each case being 
removed by means of least squares estimation 
method from the phase real power variable P, from 

15 the phase wattless component variable Q and from 

the squared phase current amplitude variable I^, 
and 

phase resistance values R being determined from 
R=P/I^ and phase reactance values X being 
20 determined from X=Q/I^, and phase impedance values 

Z=R+jX being thus determined. 

When forming the phase impedance values, it is 
particularly advantageous to remove the system 

25 frequency components (for example 5 0 Hz components) 
from the phase real power variable P, from the phase 
wattless component variable Q and from the squared 
phase current amplitude variable by means a least 
squares estimation method in each case. Such system 

30 frequency interference components would adversely 
affect the evaluation of the phase impedance values 
determined from these variables. 

In order to explain the invention further. Figure 1 
35 shows a block diagram of an exemplary embodiment of the 
method according to the invention. Figure 2 shows a 
block diagram for determining the phase impedance 
values. 



- 6 - 

Figure 3 shows the impulse response of the filters used 
for impedance determination. 

Figure 4 shows the real power and wattless component 
5 variable profiles before filtering, and 

Figure 5 shows the real power and wattless component 
variable profiles after filtering. 

10 Figure 1 shows, schematically, a method for determining 
the oscillation behavior of a three-phase electrical 
power supply system, by means of which a dedicated 
oscillation signal Pdl, Pd2 and Pd3 is produced for 
each phase in the power supply system. This is done by 
15 providing three changeover switches Ul, U2 and U3 and 
three memory elements Spi , Sp2 and Sp3 ; the connecting 
lines between the individual units in the layout are 
designed with three poles. The phase current and phase 
voltage sample values, i and u of all three phases are 
2 0 supplied to a unit for impedance determination lb, at 
whose output phase impedance values Z are output for 
the three phases. These phase impedance values Z are 
supplied via the changeover switches Ul, U2 and U3 to 
an oscillation identification unit Pe . The oscillation 
25 identification unit Pe uses the time profiles of the 
phase impedance values Z to identify the occurrence of 
any oscillation in the individual phases, for example 
in phase 1, and emits an oscillation set signal Ps at 
its output for each phase in which oscillation is 
30 identified, for example for the phase 1. The 
oscillation set signal Ps sets the memory element 
associated with the respective phase, for example Spl, 
and this memory element emits at its output the phase- 
specific oscillation signal, for example, Pdl. When an 
35 oscillation signal is being . emitted, for example the 
oscillation signal Pdl, the changeover switch, for 
example Ul, associated with the respective phase is 
switched over. The phase impedance values Z which are 
still formed for the phase in which the oscillation has 
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been identified, for example the phase 1, are supplied 
to an oscillation signal resetting unit Pu. This 
oscillation signal resetting unit Pu identifies that 
the oscillation has stopped and, in this case, emits an 
5 oscillation reset signal Pr at its output, which resets 
the memory element for the respective phase, for 
example Spl . The oscillation signal for the respective 
phase, for example Pdl, is thus also no longer emitted, 
and the respective changeover switch, for example Ul , 

10 moves back to its original position once again. In 
response to a stimulus, a unit for phase selection Pa 
ensures that the phase impedance values of the phases 
to be investigated for oscillation are in each case 
processed by the oscillation identification unit Pe and 

15 by the oscillation signal resetting unit Pu. 

The method of operation of the four units comprising 
impedance determination lb, the oscillation 

identification unit Pe, the oscillation signal 
20 resetting unit Pu and phase selection Pa will be 
explained in more detail in the following text. 

As shown in Figure 2, the phase current and phase 
voltage sample values i and u are filtered in the 

2 5 impedance determination unit lb by means of orthogonal 

FIR filters Fl, F2 , F3 and F4, thus resulting in the 
production of a variable U_re containing the real part 
of the phase voltage sample values, a variable U_im 
containing the imaginary part of the phase voltage 

3 0 sample values, a variable I_re containing the real part 

of the phase current sample values, and the variable 
I_im containing the imaginary part of the phase current 
sample values . 

3 5 Figure 3 shows the impulse responses of the filters Fl 
to F4, with the impulse response of the filters Fl and^ 
F3 which determine the real parts being annotated "o", 
and the impulse response of the filters F2 and F4 which 
determine the imaginary parts, being annotated "+". 
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As shown in Figure 2, following this, a phase real 
power variable P is calculated in accordance with 
equation (1) below and a phase wattless component 
5 variable Q is calculated in accordance with equation 
(2) in the unit 5, and a squared phase current 
amplitude variable is calculated in accordance with 
equation (3) in the unit 6. 

p = u_re-I_re - U_im-I_um (1) 
Q = U_im-I_re + U_im-I_re (2) 
= i_re ■ I_re + I_im- I_im (3) 

10 

After this, the phase real power variable P, the phase 
wattless component variable Q and the squared phase 
current amplitude variable are filtered in units 7 
and 8 in order to remove the interference 50 Hz 
15 components contained in these variables; this results 
in the filtered variables P\ Q' and The least 

squares estimation method used for this, filtering will 
be explained in detail further below. 

20 The upper illustration a) in Figure 4 shows the profile 
of the real power variable P, and the lower 
illustration b) shows the profile of the Wattless 
component variable Q, before filtering by means of the 
least squares estimation method, in each case plotted 

25 against the time t. 

The upper illustration a) in Figure 5 shows the profile 
of the real power variable P'/ and the lower 
illustration b) shows the profile of the Wattless 
30 component variable Q\ after filtering by means of the 
least squares estimation method; it can clearly be seen- 
that the 50 Hz components have been removed. 

As shown in Figure 2, after the filtering in the unit 
35 9, phase resistance values R and phase reactance values 
X are determined in accordance with equation (4) , and 
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the phase impedance values Z=R+jX .determined from them 
are emitted at the output of the impedance 
determination unit lb. 

R=P*/I^' X=Q'/I^' (4) 

5 

A least squares estimation method using a signal model 
in accordance with equation (5) is applied separately 
to each of the variables Q and in order to filter 
out the 50 Hz components contained in the phase real 
10 power variable P, in the phase Wattless component 
variable Q and in the squared phase current amplitude 
variable . 



y, = A e ' s^CDok T^j + B-e * cos{cOok T^) + C (5) 

15 

The estimation method uses the squared phase current 
amplitude variable I^, the phase real power variable P 
and the phase Wattless component variable Q to 

0 

calculate the parameters A, B and C in the signal 
20 model. The parameter C gives, the sought magnitude of 
the phase real power variable P', of the phase Wattless 
component variable Q' and of the squared phase current 
amplitude variable I^'. The summands with the 
parameters A and B model the 50 Hz components. The 
25 variable coo is the frequency (50 Hz) to be filtered 
out, and Ta is the sampling time. 

If an equivalent circuit with only two generator 
machines at the ends of a power transmission line is 

30 considered for the power supply system, the amplitude 
of the 50 Hz components decays with the time constants 
X of the sum impedance between the two generator 
machines in accordance with equation (6), where L is 
the loop inductance and R is the loop resistance of the 

35 circuit which is closed via the two generator machines. 
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(6) 



The coefficients A, B and C are determined such that 
the sum of the squares of the errors between values y 
5 determined from the phase current and phase voltage 
sample values i and u and the sample values yk 
calculated in accordance with equation (5) becomes a 
minimum (see equation (7) ) . 



In equation (7) , J represents the Q-criterion to be 
minimized. The signal model included in equation (5) is 
used as the function h{0k) • The parameters A, B and C 
15 to be determined form a vector ©k in accordance with 
equation (8) - 



20 The Q-criterion J is derived based on the parameter 
vector ©k in order, to solve the minimization task. For 
the signal model in accordance with equation (5) , this 
then results in equation (9) together with equation 
(10). 



(7) 



10 




(8) 



25 



0= i:2Y:(y,-x,©o 




(10) 
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If equation (9) is solved for the vector ©k, then this 

results in equation (11) , by means of which, and using 

the matrix Sic included in equations (12) and (13) , the 
vector ©k is determined. 



(11) 



(12) 



(13) 



Of the parameters A, B and C contained in the vector 
0jc, only the parameter C is evaluated. The vectors j:* 

10 in accordance with equation (10) and the matrix Sk in 
accordance with equation (13) are calculated and are 
stored as constants, so that they are available every 
time the method is used. 

15 Monotony criteria are applied to the locus curves' of 
the impedance values in the impedance plane in the 
oscillation identification \init Pe, in order to 
identify the oscillation process. This method for 
identification of the oscillation process is known per 

20 se, and is described in German Patent DE 197 46 719 CI. 



The oscillation signal resetting unit Pu determines 
whether an oscillation which has already been 
identified is still continuing. The procedure used for 
25 this purpose comprises the production of an oscillation 
model for phase impedance values Z associated with the 
oscillation, A check is then carried out to determine 
whether the locus curve which is described by the newly 
determined phase impedance values Z still corresponds 
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to the oscillation model. When producing the 
oscillation model, it is assumed that the locus curve 
is free of discontinuities, and that its direction 
changes only very slowly. In the present exemplary 
5 embodiment, the locus curve is described by a first 
order power function, that is to say by a linear 
equation, in accordance with equation (14) . 

X(R)=m-R+Xo ^'^^^ 

10 The parameters m and XO are determined by means of a 
non-recursive least squares estimation method from the 
last N determined phase impedance values Z. 

The linear equation is used as a model rule for the 
15 least squares estimation method, with the parameter m 
characterizing the gradient, and the parameter XO 
characterizing the offset of the linear equation. The 
parameters m and XO for the model, in accordance with 
equation (14) are determined from the last determined 
20 value pairs (Ri, Xi) of the phase impedance values Zi 
such that the sum of the squares of the errors between 
the values Xi determined from the measured phase 
current and phase voltage sample values i and u and the 
values X calculated in accordance with equation (14) is 
25 a minimum (see equation (15)). 



In equation (15) , J is the Q-criterion to be minimized, 
the model rule in accordance with equation (14) is used 
30 as the function h(0k) . In Accordance with equation 
(16) , the parameter vector ®x contains the parameters m 
and XO, to be determined, from the model rule. 



(15) 




<16) 



35 
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In order to solve the minimization task, the Q- 

criterion J must be derived on the basis of the 

parameter vector ©k- The following equations (17) and 

(18) are then obtained for the signal model in 
accordance with equation (14) : 

0= i22:rU.-A(®.X) <"> 



(18) 



If equation (17) is solved for the parameter vector ©k, 
10 then this results in equation (19) for determining the 
parameter vector 0^. 



where 



(21) 



After substitution of the parameters in equation (14) , 
this thus results in the estimated oscillation model - 
If a newly determined phase impedance value Z 
corresponds to the oscillation model, that is to say it 
20 is within a tolerance band around the linear equation 
represented by equation (14), then this identifies that 
the oscillation is continuing. If the newly determined 
phase impedance value Z is outside the tolerance band, 
then this identifies that the oscillation has stopped, 
25 and an oscillation reset signal Pr is emitted at the 
output of the oscillation signal, resetting unit Pu for 
the memory element Spl, Sp2 or Sp3 for the respective 
phase . 



30 The phase selection unit Pa receives a stimulus from, 
for example, a distance protection device which is not 
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illustrated. Depending on the nature of the stimulated 
loop, it determines the phases for which the 
oscillation identification unit Pe and/or the 
oscillation signal resetting unit Pii should investigate 
5 the oscillation behavior. The following table shows the 
association between the stimulated loops and the 
phases . 



Stimulated loops 


Phases to be investigated 
for oscillation behavior 


LIE 


LI 


L2E 


L2 


L3E 


L3 


L12 


LI and L2 


L23 


L2 and L3 


L31 


LI and L3 



